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Event by event EPOS calculations in which the expansion of the system is described by ideal 
hydrodynamics reproduce well the measured centrality dependence of V2/e par t, although it has 
been claimed that only viscous hydrodynamics can reproduce these data. This is due to the core 
- corona effect which manifests itself in the initial condition of the hydrodynamical expansion. 
The centrality dependence of V2/e part can be understood in the recently advanced core-corona 
model, a simple parameter free EPOS inspired model to describe the centrality dependence of 
different observables from SPS to RHIC energies. This model has already been successfully applied 
to understand the centrality dependence of multiplicities and of the average transverse momentum 
of identified particles. 

PACS numbers: 



Ever since Drescher et al. [T| have claimed that the elliptic flow, — , observed in the RHIC experiments, is even for 
the most central collisions at least 25% below the ideal hydrodynamical limit, the centrality dependence of the elliptic 
flow is considered as one of the key variables for the understanding of ultrarclativistic heavy ion collisions. Before this 
finding one had concluded from other observables that the Quark Gluon Plasma, created in the energetic reaction, is 
a perfect liquid @]. 

The analysis in [l[ is based on the assumption that the centrality dependence of the elliptic flow can be described 
by a simple formula 3] : 



hydro -, 
V 2 V 7 U 1 



(1) 



V2 is the measured elliptic flow, v 2 —< cos2(<f) — >, where ip (4>r) is the azimuthal angle of the emitted particle 
(reaction plane). v 2 v r ° , the value of the elliptic flow expected in the hydrodynamical limit, and Kq are the two 
free parameters of the model. The average is performed over many events which belong to the same centrality 
class of the reaction, e is the spatial participant eccentricity of the particles in an event at the beginning of the 

hydrodynamical expansion, defined as e = - a 2°^2 +4 ° x " with a 2 =< x 2 > — < x > 2 , a 2 =< y 2 > — < y > 2 , 
o~ X y =< xy > — < x >< y >. It has to be calculated in a theoretical model. K is the Knudsen number defined as 

1 a dN 

K=S^ (2) 



with a being the partonic cross section, c s being the sound velocity and S being the transverse area = 47ry a 2 a 2 — a 2 y . 

Using a Glauber model to determine S and e, the measured Q (phobos) or calculated values [|[ (star) values of dn/dn 
as well as the experimentally observed v 2 values of the Phobos Q, the Star @,[8| and the Phenix collaboration (2ll.l22| 
one can calculate ^ as a function of dn/dn. The ratio between dn/drj and dn/dy depends on the particle species but 
is not far from 1. We assume here dn/dn = dn/dy because our conclusions do not depend on the precise value of this 
ratio. The result of this analysis is shown in fig. [1] 

On the first view it seems that the data of different colliding systems at different energies fall in this representation 
on top of each other and therefore the Knudsen number depends indeed only on However, if one regards the left 

figure in more detail one realizes that at a given centrality — is almost independent of the beam energy Q. dn/dn for 
200 AGeV is about 1.4 times larger than that observed for 62 AGeV, almost independent of the centrality Thus 
the Knudsen number extracted from the 200 AGeV data is 40% smaller than that from the 62 AGeV data. This 40% 
change of the Knudsen number is hardly visible in the plot which means that this presentation is not very sensitive 
to the Knudsen number. That for the same value of dn/dn the value of — of the 62 AGeV data is systematically 
higher than that of the 200 AGeV data can be attributed to different sound velocities @, however for the price that 
an additional free parameter is employed. The speed of sound depends on the density. The density profile is different 
at different energies but approximately the same for different systems at different energies. It remains to be explained 
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FIG. 1: 



Dependence of V2 / e par t on -g ^ in AuAu and CuCu collisions at 200 AGeV in comparison with the results of the STAR 



experiment @,H and for 200 AGeV/62 AGeV in comparison with the Phobos data 



why at the lower beam energy the system is closer to the hydrodynamical limit. However, the error bars are too large 
to allow for firm conclusions. The scaling of ^ with dn/drj is questioned by the more recent Star data, displayed on 
the right hand side, where only statistical errors are shown. The Star collaboration has analyzed also peripheral AuAu 
events, in contradistinction to the Phobos collaboration. Therefore these data allow a more detailed comparison of 
both systems. We see for both systems an almost linear dependence of — with dn/drj. The slope for the CuCu system 
is, however, clearly different from that for AuAu, what one could have guessed already from the Phobos data p}. 

By fitting the experimental results with the theoretical curve, eq. [H one can determine K/(ac s ) and finally the 
Knudsen number K of each centrality bin, assuming (centrality independent) values for a and c s . Even for the most 
central bins — depends on dn/drj. Thus the hydrodynamical limit K=0 is not reached. From the extrapolation 
towards dn/drj — > oo, corresponding to K — > 0, the authors concluded that the observed flow is at least 25 % below 
the limit for ideal hydrodynamics. 

These findings triggered the development of codes using viscous hydrodynamics. Luzum et al.jioj showed in viscous 
hydrodynamical calculations that different viscosities yield a different centrality dependence of v^j 'epart and of < px >■ 
This dependence they used to identify that viscosity value which describes the 200 AGeV AuAu data. More recently, 



Song et al. [11| have extended the calculation to the CuCu system. 

What has not been realized in this line of arguments is the importance of the initial configuration at the moment 
when the hydrodynamical expansion starts. The equations of ideal hydrodynamics predict how - for a given equation 
of state - this initial configuration develops in time under the assumption that the system is locally in equilibrium. 
Therefore it is expected that different initial conditions give different final distributions. The initial configuration 
at the beginning of the hydrodynamical expansion is unknown. It is difficult to asses because the time scale for 
equilibration in the initial state is much faster than predicted by microscopical calculations. 

Recently one realized that event by event fluctuations of the initial conditions are not only much stronger than 
initially expected but also that they manifest themselves in the observables [IH, [l3[ • Only the particles deep inside 
the interaction region, the core particles, collide sufficiently frequent to form a locally equilibrated source whereas 
those close to the surface, called corona particles, do not come to a local equilibrium and do therefore not take part 
in the hydrodynamical evolution of the system. 

EPOS is a consistent quantum mechanical multiple scattering approach based on partons and strings [3], where 
cross sections and the particle production are calculated consistently, taking into account energy conservation in both 
cases (unlike other models where energy conservation is not considered for cross section calculations [HI). A special 
feature is the explicit treatment of projectile and target remnants, leading to a very good description of baryon and 
antibaryon production as measured in proton-proton collisions at 158 GeV at CERN [la ]. Nuclear effects related to 
CRONIN transverse momentum broadening, parton saturation, and screening have been introduced into EPOS [l7| . 

In heavy ion collisions (and more recently also in proton-proton collisions) collective behavior is taken into ac- 
count [18], in the following fashion: the initial scatterings, as described above, lead to the formation of strings, which 
break into segments, usually identified with hadrons. When it comes to heavy ion collisions, the procedure is modified: 
one considers the situation at an early proper time tq, long before the hadrons are formed: one distinguishes between 
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string segments in dense areas (more than some critical density po segments per unit volume), from those in low 
density areas. The high density areas are referred to as core, the low density areas as corona [lg]. It is important to 
note that initial conditions from EPOS are based on strings, providing a "flux-tube" like structure in case of individual 
events (a single flux tube in case of many overlayed events) . Based on the four- momenta of the string segments which 
constitute the core, we compute the energy density e(tq,x) and the flow velocity v(tq,x). 

Having fixed the initial conditions, the system evolves according to the equations of ideal hydrodynamics. Hadroniza- 
tion occurs finally based on the Cooper-Frye prescription. The centrality dependence of — of the EPOS calculation 
for AuAu at 200 AGeV is compared with the experimental results in fig. [5J We see an almost perfect agreement 
despite of the fact that EPOS uses ideal hydrodynamics. Especially the absolute value of v^ dro /epsilon, the free 
parameter in the core corona model, is well reproduced. The centrality dependence is due to the core-corona effect in 
the initial configuration [l8j], which manifests itself also in other observables. 



0.35 



o 
o 

ci 

> 



i 0.3 
0.25 
0.2 
0.15 
0.1 
0.05 




EPOS 2 



with hadron cascade 




w/o hadron cascade 
data: phobos 

50 100 150 200 250 300 350 400 

Npart 



FIG. 2: V2 Impart for AuAu collsions at y/s = 200GeV , as predicted by EPOS, compared with data. The curves are the result 
of event - by event hydrodynamical calculations with a freeze out temperature of 166 MeV with (solid) and without (dashed) 
hadronic rescattering. 

Recently we have developed a simple model [ill to study the consequences of this core-corona effect and to 
interpret the very involved EPOS simulations. In this model we define corona particles as those nucleons which have 
only one initial collision whereas the others are considered as core particles. / core is the fraction of core nucleons which 
depends on the centrality, the system size and (weakly) on the beam energy. This fraction is determined in the same 
Glauber model which is used to calculate the eccentricity e and the transverse area S. Core particles form a locally 
equilibrated source whereas corona nucleons are treated like elementary proton-proton collisions with no interactions 
with their environment. Of course this is a very crude model but it has the advantage that no free parameter is 
needed in this superposition of pp collisions and a thermalized source. The present experimental error bars do leave 
little room for an improvement as has been shown in fl9l . [20| . 

In this simple model we could show that, independent of the system size, the centrality dependence of the multiplicity 
M % and of the < p l T > of all identified hadrons from SPS to RHIC energies can quantitatively be described by 

M i (JV part ) (3) 

= -Apart [/core ' -^corc "F (1 — J core) ' ^corona] 

< Pt > (^part) (4) 
= Apart [/core' < Pt core > +(1 — /core)' < Pt corona > ] • 
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For our calculation, we fix M core (< p l T coro >) by applying eq. [3J Qto the most central AuAu or PbPb data point, 
-^corona an d < Pt corona > are gi ven as half 01 the multiplicity and the < p l T > measured in pp collisions. Once these 
parameters are fixed, the centrality dependence of M l and < p l T > is determined by eq. [3] and |U Especially the 
centrality dependence of the observables for the lighter CuCu system follows then without any further input. Both, the 




FIG. 3: V2 UparAN V art) / f cor e{N part) as a function of N par t for the data at 62 (top) and 200 AGeV (bottom). The data are 
fromrefs. @? ,[|| (Star) , [| (Phobos) and [UIH (Phenix). 

core - corona model as well as the approach of Drescher et al. pj , assume that the system is not in the hydrodynamical 
limit, however the reason is different. In the approach of [l[ it is assumed that all particles come to the same degree of 
thermalization independent of whether they are at the surface or in the interior of the interaction zone and that this 
thermalization is not complete. The degree of thermalization depends on the centrality. In the core - corona model 
the particles in the interior of the reaction zone have sufficient collisions to reach complete thermalization whereas 
those located at the surface do not came to equilibrium at all. They behave like being created in independent NN 
collisions. In this separation it follows the findings of the EPOS program. The centrality dependence is due to the 
centrality dependence of f core- 
In the core - corona model only the core particles feel the eccentricity of the overlap region. The particles produced 
by corona nucleons show an isotropic distribution in azimuthal direction. Therefore, in the core - corona model, 
V2 Impart is expected to be oc f CO re 

hydro 

1>2/ ' £part{N par t) = fcore(Npart)-^—^ ■ (5) 

Because f C ore(N par t) and e par t is calculated in the Glauber model v 2 a r ° is the only free parameter. The core - 
corona model allows for some immediate predictions. Because f cor e is similar in central CuCu and AuAu collisions 
(Fig. [3]) we expect that central collisions of AuAu and CuCu show a similar v-ijtpart despite of the large difference of 
Npart ■ Therefore V2 / e par t plotted as a function of ^ ^ is in central CuCu data larger as compared to AuAu data. In 
peripheral CuCu and AuAu collision V2 / (-part should be very similar for the same N par t ■ The easiest way to compare 
our results with experiments is to divide the experimental results by f C ore{N pa rt)- In this presentation we divided out 
the centrality dependence and therefore we expect that the data arc horizontal straight lines. 

In Fig. [3] all presently available experimental results are presented in that way. We see that in this presentation 
the 200 AGeV data of Phenix and Phobos, presented in the bottom part, follow a straight line. The Star data show 
a small slope. We should stress that for the peripheral CuCu data where the core contains about 15 nucleons we 
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cannot expect that our assumptions are fully justified and one may doubt whether there the core forms really a 
fully thermalized source. Also the data at 62 AGeV, shown in the top panels, are compatible with our prediction 
that V2 Vdro /e P art = V2 / £part(N pa rt) / fcore is constant. However, whereas the Phobos data are compatible with the 
assumptions that V2 Vdro /e part is the same for AuAu and CuCu, the Phenix data show a difference of 30% between 
CuCu and AuAu. This is not understood yet. 

Our model can be extended to the lower SPS energies where V2 / e par t(Np ar i) has been measured by the NA49 
collaboration (23[. These results, show in fig. [H are also compatible with a straight line. 
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FIG. 4: V2/ e P art (Npart)/ fcore(N pa rt) as a function of Npart for the NA49 data [23] at y/~s = Y1.2GeV . 

In conclusions, we have shown that the experimentally observed centrality dependence of vi / ' e par t(N part ) is well 
described by the EPOS event generator using ideal hydrodynamics. In contradistinction to other hydrodynamical 
calculations EPOS creates its initial condition which fluctuates from event by event. Only those regions of the 
interaction zone which have a density above a threshold density form a locally equilibrated source whereas the others 
are treated as independent NN collisions. The simple EPOS inspired core - corona model, which has successfully 
described the centrality dependence of the multiplicity and the average transverse momentum of identified particles, 
can also quantitatively explain the centrality dependence of vij ' e pa rt(N part ) . Only core particles come to a local 
thermal equilibrium and develop a whereas the corona particles decay isotropically. The fraction of both classes of 
particles depends on the centrality and the centrality dependence of the elliptic flow reflects directly this dependence. 
The only free parameter is V2 Vdro . This results questions the idea to use the centrality dependence of V2 to determine 
the viscosity of a quark gluon plasma. First of all, on the conceptual level, it is difficult to imagine that particles 
close the surface of the interaction zone behave exactly as those in in the center. Second, on the quantitative level, 
two additional free parameters, Kq and c s (E), have to be introduced to relate the centrality dependence of i>2 with 
the expectations for a viscous system. The numerical value of both is difficult to asses by other observables. 

All approaches agree on the fact that the system is not a perfect liquid, even in central collisions. In the core corona 
model this is due to the fact that nucleons at the surface of the interaction zone behave like pp collisions whereas 
the center comes to complete equilibrium. In the viscous hydrodynamical approach it is assumed that the interaction 
zone has no surface but that the whole system is described by a finite viscosity which varies with centrality. 
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